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Abstract 

We make a phenomenological model of optical two-beam interaction in a model planar liquid 
crystal cell. The liquid crystal is subject to homeotropic anchoring at the cell walls, is surrounded 
by thin photosensitive layers, and is subject to a variable potential across the cell. These sys- 
tems are often known as liquid crystal photorefractive systems. The interference between the two 
obliquely incident beams causes a time-independent periodic modulation in electric field intensity 
in the direction transverse to the cell normal. Our model includes this field phenomenologically 
by affecting the potential at the walls of the cell. The transverse periodic surface potential causes 
spatially periodic departures from a pure homeotropic texture. The texture modulation acts as a 
grating for the incident light. The incident light is both directly transmitted and also subject to 
diffraction. The first diffracted beams result in energy exchange between the beams. We find that 
the degree of energy exchange can be strongly sensitive to the mean angle of incidence, the angle 
between the beams, and the imposed potential across the cell. We use the model to speculate about 
what factors optimize non- linear optical interaction in liquid crystalline photorefractive systems. 

PACS numbers: 42. 15.-1, 42.25.Fx, 42.70.Df, 42.70.Nq, 42.79.Dj, 42.79.Kr, 61.30.-v 
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I. INTRODUCTION 

It has long been known that nematic hquid crystals act as non-linear optical media 
The electric fields in a strong light beam reorient the liquid crystal director. In so doing they 
affect the dielectric properties of the medium, and hence its light transmission and refiection. 
Thus the liquid crystal reacts differently to a high intensity beam than when the intensity 
is low. A slab of liquid crystal exhibits analogous properties when irradiated by two beams 
rather than by a single beam. The liquid crystal responds to the interference pattern between 
the beams. The result is a grating in the liquid crystal cell, which diffracts the incoming 
beams. The lowest order diffracted beams from each incident beam act to reinforce the 
other, leading to the phenomenon of beam amplification. This latter phenomenon is the 
subject of this paper. 

In general beam-coupling is a highly non-linear phenomenon and as such requires intense 
beams to manifest itself. However, it has been found experimentally that there are circum- 
stances when beam-coupling appears to occur at much lower light intensities. The device 
possibilities of these high beam-coupling conditions have meant that these systems have 
attracted much interest. Two particular interesting systems obtain when either the liquid 
crystal is doped by dye molecules, or when the liquid crystal cell is sandwiched between 
walls consisting of photoconducting material. Now free charges can play an important role 
in the non-linear optics - in the former case the ions move inside the liquid crystal itself, 
whereas in the latter case they affect the boundary conditions to which the liquid crystal 
is subject. In this paper we consider the second of these cases - that in which the liquid 
crystal is surrounded by photosensitive layers. An associated feature of such systems is that 
the degree of beam-coupling is strongly dependent on, and amplified by, a low-frequency 
voltage across the liquid crystal cell. 

The existence of ions in these systems has led to this phenomenon being linked with 
photorefraction j^. In photorefractive systems moving charges lead to non-linear optical 
effects. These so-called photorefractive liquid crystal systems exhibit large optical non- 
linearities, at least partly because two non-linear optical processes seem to be manifesting 
themselves simultaneously. This statement, however, while true, is insufficient even to give 
the most basic description of the physics of beam-coupling in these systems. In this paper 
we present a phenomenological model, which although by no means a complete description 
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of the system, provides a basic understanding of some of the most striking features of the 
experiments. The most important of these features is the observation that the largest beam- 
couphng effects occur when the grating period is comparable to the cell thickness. 

Experimental work in photorefraction in liquid crystals dates back about decade. In the 
cases of interest in this paper, the effect results because a spatially modulated light fiek 
causes a modulation of the electric field either in the aligning layer itself 4|, js], [g], [t] 
or in the interface between the LC and the aligning layer [s], M, • In the first case the 
liquid crystal cell is lined by photoconductive aligning layers, whose electrical resistance is 
decreased by light irradiation. This increases the electric field in the liquid crystal bulk, 
which in turn causes a spatially modulated reorientation of the director in the cell. The 
effect is reversible; the induced gratings disappear when the incident light is switched off. By 
contrast, in the second case the photorefraction is controlled by the processes in the interface 
between LC and aligning surfaces. Both of these layers may be nominally insensitive to light. 
The resultant spatially modulated electric field induces a reorientation of the director in the 
bulk and a permanent grating. 

Theoretical work on these systems has concentrated on extending existing photorefractive 
concepts, which have been developed for optically isotropic systems. In this case, the key 
inputs into an experiment are: the cell thickness L, the light wave-length A, the grating 



period A, and the dielectric constant e of the isotropic medium. Kogelnik 



ll| developed 



a coupled-wave theory which can predict the response of volume holograms (i.e. thick 
gratings). Klein gave the criterion for a grating to be thick in terms of the parameter 
Q = The coupled wave theory begins to give good results when Q > 10. 



Montemezzani and Zgonik [13| have extended the Kogelnik coupled-wave theory to the 
case of moderately absorbing thick anisotropic materials with grating vector and medium 
boundaries arbitrary oriented with respect to the main axes of the optical indicatrix. The 
dielectric tensor modulation takes the form 

e= [e° el cos {K ■ r)] +i [5° + e] cos {K ■ r + (p)] . (1) 

Galstyan et al jl^ have also presented a variant of this idea applied to anisotropic thin 
holographic media. 

The key extra piece of physics in liquid crystal cells is that the director is anchored by 
the cell walls. As a result the spatial modulation of the dielectric function is considerably 
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more complicated than the Montemezzani-Zgonik form. In addition, the hquid crystal cell 
parameters are often in the so-called Raman-Nath regime Is], [l^, which corresponds to thin 
gratings. For thin isotropic gratings with a one-dimensional refractive index modulation, 



;he theory is well developed (see for example [15[). For such a system, for example, Kojima 
~v\ used a phase function method to understand the diffraction problem for weakly inho- 
mogeneous anisotropic materials in the Raman-Nath regime, assuming a dielectric function 
spatial modulation e ~ cos {Qt — Kx). 

In this paper we shall study the diffraction and energy transfer of two light beams inter- 
secting in a nematic liquid crystal cell with strong homeotropic anchoring at the cell walls. 
The beam-coupling can be amplified by a DC-electric field, which is applied to the cell 
perpendicular to the cell walls ((92;-direction). This problem corresponds to that addressed 
experimentally by Korneichuk et al 1^. 

The presence of the two beams causes a periodic lattice in the light intensity field in the 
cell bulk and its boundaries. In addition, the laterally periodic light intensity also causes 
a modulation in the dc-electric field potential at the cell boundaries. This paper addresses 
the photorefraction problem phenomenologically. 

Specifically, we consider here a restricted problem with two major caveats. Firstly we 
shall not examine too closely the origin and mechanism of this modulation. We simply 
remark that it can and does result from different physico-chemical phenomena taking place 
at the cell walls. Likewise in this simple approach, we shall suppose that the physics of the 
electric field in the liquid crystal is driven by dielectric processes, and that charge transport 
does not play a major role in determining director orientation or light scattering. Elsewhere 
we shall relax both of these constraints. 

The key to understanding beam-coupling in these systems lies in the following observation. 
The surface potential modulation produces a spatially modulated electric field. The resulting 
torque on the liquid crystal director distorts the initial homogeneous homeotropic alignment. 
The consequence is an anisotropic medium with a spatially modulated director and hence 
optical axis. The test beam - or the beams that write the grating - diffract from the liquid 
crystal cell, which now possesses a spatially modulated refractive index. One may then 
calculate beam diffraction and inter-beam energy transfer. 

The paper is organized as follows. In Section II we determine the electric field profile in 
a cell subject to a light-induced periodic modulation of the surface potential. In Section III 
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we calculate the director distribution inside the liquid crystal cell subject to this spatially 
modulated electric field. Then in Section IV we present results of calculations of beam 
diffraction and energy transfer. Finally in Section V we present some brief conclusions, and 
focus on possible extensions of the model. 




FIG. 1: Schematic picture of a two-beam coupling experiment, showing meaning of quantities used 
in the paper. 



II. ELECTRIC FIELD WITHIN THE 
LIQUID CRYSTAL SLAB 

We consider two equal frequency light beams with wave numbers ki ^ k2 inside the 
medium. These beams give rise to electric fields Ei, E2, and intensities /i, I2 proportional 
to the squares of the respective electric fields: 



El = Eioexp(zki -r), (2) 

E2 = E2oexp(2k2-r), (3) 

ki = (/csin + 7) ,0, fccos (t/^ + 7)) , (4) 

k2 = (/c sin (z/; — 7) , 0, cos (z/^ — 7)) . (5) 

The bisector of the beams makes an angle ^^ with the cell normal and 2'y << ip is the 

angle between the beams. Initially we suppose that the scattering by the cell is weak, and 
thus the light transmission through the cell is close to unity. 
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The beams interfere in the hquid crystal slab, forming a complex intensity pattern with 
wave number q = ki^ — k2x- In principle the pattern must be calculates self-consistently. 
In practice it may be possible to treat surface-induced and bulk-induced effects separately. 
In this section we concentrate on the particular effect at the surfaces. At the bottom (i.e. 
incident) interface, the interference pattern of light takes the form: 

I{x,z = 0) = Ji + /2 + 2^(/i/2)i/2cosga;, 

q = kix- k2x (6) 

Likewise, at the top substrate we have an analogous pattern, but shifted in phase with 
respect to the lower substrate: 

I{x,z = L) = h + I2 + 2 Y.{hhY'^ cos{qx + 5), 

q = kicc- k2x (7) 

6={k,,-k2,)L (8) 

In the absence of the light beams, we suppose a voltage $0 across the liquid crystal cell. 
This is a key input to the theory. In photorefractive systems the optical non-linear effects 
are large and strongly amplified by a voltage across the cell. In the simple theory presented 
here the non-linear effect in the absence of an external field is strictly zero. 

We now make the hypothesis that the spatial distribution of light intensity induces a 
modulation in the surface potentials. The effect of the light beams is to modify these 
potentials slightly. The boundary conditions on the electric potential at the top and bottom 
substrates can now be written: 



(p{z = 0,x) = -$0/2 + $1 cosgx, (9) 
ip{z = L,x) = $o/2 + $icos(ga;-h5) (10) 

where $1 = a(/i/2)^''^. The parameter a is a phenomenological quantity, and in principle is 
different for each surface configuration. Here we are supposing that the surface preparation 
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of the upper and lower surfaces is identical. 



independently by a Frederiks experiment 
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'he parameter a can in principle be determined 



We now proceed to determine the electric field potential within the liquid crystal slab. 
The electric field obeys the equation 

V-D = (11) 

with Di = SijEj = [e±Sij + (en — n^nj] Ej, where n is the nematic director. 

We can solve eq. ffTT]) using relation E = — Vy?. At this stage we note that the equations for 
ip and n must be solved self-consistently. However the liquid crystal is subject to homeotropic 
boundary conditions, and hence except at very strong light intensities the director is closely 
aligned to the direction perpendicular to the slab: n ^ e^, and then 



The problem to be solved is thus eq.( ]T2i) . subject to the boundary conditions eq. ljlO 
This problem can be solved analytically, yielding 



, , ^ z-L/2 ^ 



cos 5 — cosh (gL) . , ^1 sm6 

cosh [qzj-\ — — smh [qzj> cos qx — — smh [qzj sm qx 

smh qL j smh qL 



(13) 

where q = . —q- The potential in the liquid crystal slab consists of the externally imposed 
voltage plus a contribution linear in the surface perturbation induced by the light-beam 
interference. This perturbation has the same periodicity in the direction in the cell plane 
as the initial perturbation. However, the behavior is complex as a result of the competing 
effects of the out-of-phase surface perturbations. 

The resulting electric field, in addition to the imposed external field £"0 normal to the 
cell, has components in both the x and z directions. The contributions in the x direction are 
particularly important, because they lead to director distortion, and thus to refractive index 
modulation. It will be this refractive index modulation which induces the beam-coupling 
which we seek to describe. 

The electric field inside the cell bulk is given by taking the gradient of eq. (fT3!) . We find: 
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Ex (x, z) 
Ez {x, z) 
Ei fx, z) 



Eiq {z) COS qx + E2q (z) sin qx 
Eo + K 

E^q{z) COS qx + E^q {z) sin qx 



(14) 



where 



$0 

-C/n — ^, 



(15) 



sin 5 _ 
smn qL 

„ ^ 1 ~ cos 5 — cosh qL ^ 
E2q = q'Pi I coshgz H . ^ „^ smhg^; 



3q 



sinh gL 



cos 6 — cosh gL „ 

-gOi smh H — — cosh qz 

sinh qL 



(16) 
(17) 
(18) 



4g 



sin 5 



cosh qz 



(19) 



sinh gL 

It will be useful later to normalize the field Eiq with respect to the externally imposed field: 



eiq [Z) 



Ei ^ ^ 

— ^ = Oj cosh qz + bi sinh qz, 
Eq 





fej are given by: 




i 




h 


1 





^ $1 sin 5 
$0 sinh qL 


2 


^^^^ 


$1 cos (5 — cosh qL 

n J 

$0 sinh qL 


3 


_^ $1 cos 5 — cosh qL 
$0 sinh gL 


$0 


4 


sin 5 
$0 sinh qL 






(20) 
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III. DIRECTOR PROFILE IN THE LIQUID CRYSTAL CELL 



We now determine the director profile in the hquid crystal cell in the presence of the 
electric fields given by eqs. fll4l) . The bulk free energy Fy of a distorted nematic liquid 
crystal in an applied electric field takes the form: 

Fv = ^Ku J (V ■ n)' dV+^K22 J (n ■ V x n)' dV+^K^s J (n x V x n)' dV-^ J D ■ EdV 

(21) 

with the electric displacement D = eE, Sij = e±6ij + eaUinj, with the anisotropic part of 
the static dielectric constant Sa = £\\ — £±- 

The electric field felt by the liquid crystal molecules has a number of contributions. The 
first is the externally-imposed voltage. The second is the periodic modulation in the x 
direction discussed in the last section. This is an indirect effect of the light field acting on 
the surface layer, transmitted into the bulk as a result of the effect of the Laplace equation. 
A final contribution comes from the direct effect of the light field on the liquid crystal. We 
are assuming here that this can be neglected. The justification for this is empirical, and 
derives from the observation that in the absence of the surface layers, the effect essentially 
disappears 0, isl. The director field is now given by n = (sin 6' {x, z) , 0, cos 6* (x, z)), with 9 
small. 

The variational problem to be solved consists of minimizing eq. fl2ip . subject to strong 
anchoring homeotropic boundary conditions 6 [x, z = 0) = 6 {x, z = L) = at each wall, and 
subject also to the electric fields given in eq. (fT4|) . We simplify further by supposing the so- 
called one constant approximation, i.e. the splay and bend Frank-Oseen elastic coefficients 
are equal: Kn = K33 = K. 

The relevant part of the thermodynamic functional eq.( l2T]) is now given by: 



ea{{El - E^) sin^ 9 + E^E, sin 29)]dxdz. (22) 



The Euler-Lagrange equation for this functional is: 



Sa [{El - El) sin 9 cos 9 + E^E^ cos 29) = 0. (23) 
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E E' 

We solve eg. (1251) in the limit — ^ < 1, ^ < 1- This corresponds, roughly speaking, to 

Eq Eq 

Ex 

high voltage or low beam intensities. Expanding to linear order in 9 and — , we obtain: 

Eq 



+ 



(9x^ dz 



e 



Ex 

En 



(24) 



where the length scale ^ is the relaxation length set by the bulk electric field, with ^ = 

e E"^ 

■ In order to do this, we linearize eq. fl24p for components 6*1^ = ^i{z) cosqx and 

02q = "^2 (z) singx of the director reorientation of wave number q in the cell plane. 

It is convenient at this stage to reformulate the problem in terms of non-dimensional 

variables. We define a rescaled length along the cell p = — , p G [0, 1], a rescaled transverse 

^ ij 

wave vector /i = gL, and a rescaled voltage z/ = — = LEq^-^Y^"^ . For some purposes it 
is convenient to measure length not from the plane of incidence, but rather from the mid- 
plane of the cell. We can define a length variable a = p — 1/2, and then inside the cell 
a e [-1/2,1/2]. 

Eq. flM|) now reduces to 



with = /i^ + 



dH, (p) 

dp2 



- K^di (p) = -v'^Ciq (p) 



(25) 



This equation can be solved using standard methods and has solution: 



$1 

9{x, z) = —qL— ( cos qx sin 5 

$0 



sinh pp sinh^p 
sinh p sinh k, 
cos 6 — cosh p 



+ sin qx 



cosh pp — cosh Kp + (26) 



. _ cosh K — cos 5 . , 
smh pp H — smh up 



sinh p sinh k 

For some purposes it is more convenient to rewrite eq. (l271) as: 



e(x,z) 



-qL 



$1 



COS 2 sm 



sm cos 



{qx + I) 



cosh /jcr 
cosh 



cosh Kg 
cosh ft/2 



+ 

(riT -I- -\ j sinh /JO- _ sinh /-co- \ 
yqx 2) ]_sinh/j/2 sinhK/2J 



(27) 



2 Vi"" ' 2/1 sinh/j/2 

The advantage of the expression (1271) is that the variation of 9{x, z) is expressed in terms 
of components that are respectively out-of-phase and in-phase with the total optical field 
intensities on the mid-plane. This expression explicitly exhibits the symmetry of the system 
around the mid-plane. 
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We discuss in the next section in detail how to use 6{x, z) to calculate non-linear optical 
effects. However, we note that in general the larger the values of 6{x, z) measured in some 
sense the larger will be the optical effect. It is therefore of some interest to monitor the 
behavior of 9{x, z) as a function of system parameters. 

We plot the z-dependence of the out-of-phase component of d{x, z) in FigO One might 
expect a roughly sinusoidal dependence, with a maximum at the cell mid-plane. And indeed, 
for closely matched incident beams, with a low wave-number interference pattern, this is 
what occurs. But when the non-dimensional grating wave-vector /i is larger than unity, the 
sinusoidal dependence no longer holds. By /i = 4, the response is flattened, and by ~ 6 
the profile has developed a double hump structure. Not only is the shape unexpected, but 
the magnitude is reduced in this regime, and as discussed in the last paragraph, this should 
(and, as we shall see below, does) lead to a reduced non-linear optical effects for larger /i. 

0.1 1 • ^=^^ 



0.08 - 



^0.06- 



0.04 - 



0.02 



1 
■4 

7 

■ 10 




0.5 



FIG. 2: Functional dependence of angular distortion from homeotropic texture, in the limit of low 
field, for different values of the non-dimensionalized grating wave number /i. This figure shows the 
component of the angular distortion out-of-phase with the intensity modulations. 

f I I \ 

From eq. (1271) we observe that 6 = — — — -— m~rT can be regarded a figure 

\cosh(/i/2) cosh(/c/2) / 

of merit for the degree of distortion of the liquid crystal. This is a measure of the amplitude 
of the response at the mid-plane of the cell. We note that 6 technically only measures the 
out-of-phase distortion, and furthermore even then the magnitude of the distortion is not 
maximal at the mid-plane of the cell. Nevertheless it serves in a rough and ready way as a 
surrogate for the magnitude of the grating response to optical probes. In Fig. [3] we plot G 
function of voltage. 
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FIG. 3: Dependence of the amplitude of the component of the angular distortion out-of-phase with 
the optical intensity, as a function of non-dimensionalized voltage 

IV. DIFFRACTION OF LIGHT BEAMS 
A. Formulation of problem 

We now consider light beam propagation of each of the two waves through the (now) 
weakly non-uniform anisotropic liquid crystal cell. We suppose the wave incident from the 
vacuum to have wave number k = /ck, with k = 27r/A = uj/c, where the angular frequency 
u! and the speed of light c take their usual meanings. In the presence of a uniform liquid 
crystal, the light will be refracted into an ordinary (o-) and an extraordinary (e-) wave. 
The ordinary wave is polarized perpendicularly to the plane of incidence (in the y direction) 
and the refractive index which corresponds to this wave is Uo = ^/s^- For the extraordinary 
wave the effective refractive index is given by: 



where '0 is the angle between the director and the direction of propagation inside the medium. 

The effect of the non-uniformity will be to modulate the amplitude of the wave in the 
plane of the outgoing surface, and hence in addition to the refraction, diffraction will also 
occur. The purpose of this section is to calculate the magnitude of this diffraction. There are 
in fact two waves, but first we discuss the effect of the modified director on each individual 
wave. 




(28) 
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B. The dielectric function 



The dielectric function is given by 



■J' 



(29) 



with Ea = B\\ — e± and director components rij. In the hmit of interest in this paper, the 
director deviations from the initial homeotropic alignment are small. Then 



n = [sin 9 (x, z) , 0, cos 9 {x, z)] ~ [9 {x, z) , 0, 1] 



The dielectric function now simplifies to 



e = £0 + z)ei, 



(30) 



(31) 



or alternatively: 



^ Q ^ 

ex 

£11 



+ 9{x,z) 



Q Sa^ 





(32) 



C. Geometrical Optics 

The theoretical strategy involves determining perturbations around the transmission 
through the pure homeotropic (i.e. Eq) system. The characteristic length for director inho- 
mogeneity in the ^-direction is the cell thickness L. In the x-direction the corresponding 
characteristic length is the gra ting period A = 27r/g. We shall use the Geometrical Optics 
Approximation (GOA) 
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22], valid in the limits X << L and A << A. 



We seek solutions to the Maxwell equations 

V X E = 

in the following forms: 



dH ^ d 



(33) 



E(r,t) = exp{-iujt + ikS{r))Eo; 
H(r,t) = exp{—iut + ikS{r))'H.Q 



(34) 



The term S{r) is the optical path length or eikonal, and the local direction of the wave 
vector is given by VS'(r). 
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Substituting eqs. (lMI) into the Maxwell equations (133|) . we obtain the following pair of 
equations: 

kVS{r) X Eo = /iot^Ho; kVS{r) x Hq = -eoioo'Eo. (35) 
The magnetic field Hq can now be eliminated, yielding a homogeneous equation for Eq: 



V5(r) X (V5(r) X Eq) + eEq = 0. 



(36) 



Eq. fl36p is a homogeneous system of linear equations for the electric field components, 
analogous to a vector Helmholtz equation. In general solutions to this equation will be 
trivial and uninteresting. However, there are non-trivial solutions, corresponding to optical 
traveling waves, if the determinant ofling waves, if the determinant of this set of equations 
is null. 

In fact the determinant factorizes. An E eigenvector in the y direction corresponds to 
the ordinary wave. The perturbations in the dielectric tensor do not affect transmission 
of the ordinary wave through the sample, and we shall not be interested in this mode of 
transmission. The E eigenvector in the x — z plane (i.e. the plane of incidence) corresponds 
to the extraordinary (e-) wave. The pair of homogeneous equations are: 



d,Sd,S + eJ{x,z) eu-{d.,Sf 



E^ 



E. 



(37) 



The null determinant condition appropriate to the e- wave now reduces to 29|: 



{e± - (d^Sy) (£|| - (d^^S)^) - {d^Sd^S + SaOix, z)f = 0. 



(38) 



D. Perturbation Theory 

In the absence of the director modulation, the e-wave is directly transmitted. We consider 
eq. fl38|) as a perturbation of this process. We therefore recast eq. fl38l) to lowest order in 

e{x,zy. 

{d^Sr ^ {d^ ^ 2 ( ^) ^d.s d.S) e{x, z) = 1 (39) 



In the spirit of the WKB approximation, the solution of eq. (!39l) can be expressed as the 
sum of an unperturbed e-wave, plus a small phase change /(r) which can be ascribed solely 
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to the modulation. Thus: 

5(r) = 5o(r)+/(r), (40) 

where 5*0 obeys the equation 

jd.Sof ^ jd.Sof ^ ^ ^^^^ 

We note also that the effective refractive index Ueff can be defined as follows: 

{VSoY = n%,. (42) 

Combining eqs. fl41l) and fl42p yields the well-known expression for the refractive index 0281) . 
The wave vector inside the medium is now given by 

k' = WS'o = k[d^SQ, 0, d,So] = k'[smij, 0, cos^fj], (43) 

with k' = krieff. 

Combining eqs. fl39|) . (H0!) and pTl) . we can obtain the leading order equation for the phase 
f{x,z): 

- {d.So){d.So)'-^^^^. (44) 

We now substitute eq. (H3l) into eq.(jllj), yielding: 

'sin'?/'\ df ( cost/^A d] 



dx \ £± J dz 

— sin ijj cos ip ^ ^"^"^-^-^ ^ 6 {x,z) . (45) 

Eq. (145|) can be solved using the method of characteristics. The left hand side of this 
equation can be transformed into a total derivative: 

dz dz \dz J \dx 

dz ~'~ £11 cos?/'} i^dx) dz~^^^ ^ i^dx} ^'^^^ 



Combining eqs. fl45l) and (146|) yields: 



^ = - sin (^j-^ neffOix, z), (47) 
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where now the quantities x and z in this equation are exphcitly related by 

^ = tan?/^', (48) 
dz 

defined in eq. (l46|l . Eq. (l48!l allows a family of solutions 

x{z,xq) = xq + ztaji%lj' . (49) 

We note that each member of this family of solutions represents a wave entering the liquid 
crystal sample at position xq. The direction of wave propagation is given by the angle 
However, because this medium is anisotropic, the angle of the energy propagation, given 
by the Poynting Vector, is determined by the angle ip' . The family of solutions x{z,xq) 
corresponds to paths in the undistorted anisotropic medium with different Xq travelling in 
the direction of the Poynting Vector. 

Now we can solve eq. P7|) directly, by integrating the right hand side, yielding: 

f{x,z) = -simP f£!^) r 9{x'{z',Xo),z')dz', (50) 



. ^11 / ^0 

where the integration path is such that x = Xq + ztanip'; x' = Xq + z'tan'0', and hence 

x' = X — [z — z') ta.nip'. (51) 

Thus 

f{x,z) = — sinip ^"'^^■f-f 

x[ e{x - {z - z')ta.nilj',z')dz', (52) 



The key quantity of interest is the phase retardation of the beam as it leaves the cell, i.e. 
at z = L. We now rewrite eq.( l50i) . so as to express this quantity directly: 



f{x,L) = -sinV^I^^ 

L 

X j 9{x-{L-z') tan i:', z') dz' (53) 



E. Diffraction Pattern 



The formula (153|) applies to all incident light beams. We confine our interest to the cases 
in which there are two incident beams, with wave numbers ki and /c2, with ki^ — k2x = Q- 
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Eg. (153!) permits the calculation the light fields of the beams as they exit the liquid crystal 
cell. By suitably decomposing these light field into Fourier components, it is possible to 
identify the amplitudes of particular diffracted beams. 

The light field along the plane z = L is modulated by the factor 



exp[ikS{x, L)] = exp[i(k' ■ r + A;/(r))] = 
exp[i{kxX + k'^L + kf{x, L))] = exp[{i{kxX + 50o + ^(pi)], 

where (5(/)o = k'^L = kuef jL cos ip, and 6(j)i = kf{x,L); 

6<Po = {kL cos 4j) 

COS"' + £j_ sm ipj 

and 



(54) 



(55) 



X 



v/gl(g|| - g±) sinV; 
^(^ll cos^ 'ip + £± sin^ ■i/;)i/2 

9 [x — {L — z') tan ■?/;', z) dz 



(56) 



The quantity 5(f)i{x) is the varying component of the additional phase of the incident beam, 
following from the director modulation in the liquid crystal medium. We now recall that 
from eq.( l27|) the director profile can be written in a sinusoidal form: ^(x, z) oc cos(gx + A). 
Combining this result with eq. fl56|) yields the result: 



!)i(x) = B cos{qx + A), 



(57) 



where the phase modulation parameters B > and A are respectively the amplitude and 
phase of the additional phase 50i(x): 



B = VA^ + C2; tanA = -^. 



(58) 



After complicated but straightforward algebra using eqs. (j27|) . (|56|1 . expressions for the 
quantities A and C can be derived: 



A 



-kL 



$1 fi{ei\ — e±) simp 
$0 cos^ i> + £± sin^ tjjy/'^ 



1/2 



1/2 



COS 



5 J cosh ficr cosh Kg 



2 1 cosh/i/2 coshK/2 



COS 



(/i(cr — |) tamp') 



sin f - i^} sin (/.(a - 1) tan^O 

(59) 



da 
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C = -kL 



1/2 



$1 fj,{e\\ — e±) sinip 
$0 i^w cos^ ip + £± sin^ ipy/'^ 7-1/2 



cos 



2 1 cosh /i/2 



sin I <^ - 4isi^ ^ cos 

2 1 smh/i/2 smhK/2 



(/i((T — |) tan '?/'') 

(60) 

Now eq. flMl) can be rewritten, using eq. fl57|) . so as explicitly identify different components 
of the diffracted wave. The key relation is the Jacobi-Anger expansion 23l |: 



da 



exp (iz cos 0) = (i)" J„ (z) exp 



(61) 



Combining eg. (1611) with eqs. (lM|) . (155|) . fl56|) and fl57|) yields the expansion for the electric 
field at the output surface: 



Eout = Eoexp{iS(l)o + ik^x) exp(i(50i) = 

oo 

Eq exp(i(50o + ikxx)y~^^ {i)^Jn{B) exp(inqx + inA) 



(62) 



Now we can identify 2J| terms in this expansion with the amplitudes X*^"-' and phases (5^"^ 
of outgoing waves in the diffraction pattern: 



and 



= + nA. 

The electric field in the diffracted wave of order n then takes the form: 

§^ = XW exp2(k(")-r + 5(")), 

Eq 

where k^"^ is the wave number of the diffracted wave of order n, with /ci"^ 
|k(")| = Ikl. 



(63) 
(64) 

(65) 

k^ + nq, and 



F. Beam Coupling 

We now return to the original problem ([5]) in which there are two incident waves with 
wave numbers ki, k2, and with ki^ — k2x = Q- Beam coupling corresponds the diffraction of 
waves from incident wave ki to outgoing wave k2 = ki — qe^, and from incident wave k2 to 
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outgoing wave ki = k2 + ge^.. Thus the diffracted wave of order —1 from ki adds coherently 
with the directly transmitted wave k2, and the diffracted wave of order +1 from k2 adds 
coherently with the directly transmitted wave ki. Equivalently, using the notation of the 
last section, ki = 'ki^^^ and k2 = k^"^"*. 

We are thus able to use terms from the diffraction expression eg. (1^51) to evaluate the 
amplitudes of the outgoing waves directions ki and k2. We find: 

(66) 

Eoutik2) = EoiX^-'^ exp(z5(-^)) + Eo2X(°) exp{i5^'^) 

We note that in principle the quantities X*-"-* depend on the value of the incident wave 
number. However, in the two-beam coupling case discussed here, the incident waves have 
wave numbers very nearly equal to each other, and so we may consider the quantities X^"'^ 
to be the same for each incident wave. 

From eq. fl66l) . we can evaluate the outgoing wave intensities, using the relation / = E ■ E*. 
Using eq. (p^ . We find: 

= (Eoiexp(2<5(°))Jo(5) +^Eo2exp(^5(+i))Ji(5)) 

{{E*, exp(-z5W) Jo(i?) - tE*, exp(-^5(+i)) Ji(i?))), (67) 



or 



Iout{ki)=hJo{B) + hJl{B) - 2^hhUB)J,{B) sin(A); 



(68) 



IoutiM2)=hJl{B) + hJl{B) + 2^hhUB)J^{B) sin(A), 



where A is as defined in eq.f lSSp . 

It is usual to consider one of the beams as the pump beam and the other as the probe 
beam. Without loss of generality, we shall suppose that ki corresponds to the pump beam 
and k2 to the probe beam. In line with the literature, we define 

m = ^ = (69) 

We can now rewrite the formulas for the outgoing beam intensities in terms of the quantity 
m as follows: 



20 



Iout{ki)=h (jo{B) + mJl{B) - 2^Jo{B).h{B) sin(A) 



Iout{\^2)=h [mJl{B) + Jl{B) + 2V^UB)J^{B) sin(A)] 



(70) 



The degree of beam coupling can now be characterized by the Gain g. This is the ratio 
of the intensity of the outgoing beam in the direction of the probe beam in the presence of 
the pump beam to the intensity of the same beam in the absence of the pump beam. In the 
context of this paper, in which we do not consider reflection and refraction at the cell walls, 
the quantity g is defined as: 



9 = = + -Jt{B) + 2-=MB)Jr{B) sin(A). (71) 



A related quantity is the Diffraction Efficiency rj, which measures the strength with which 
the grating diffracts the probe beam. The formal definition is the ratio of the intensity of 
the diffracted probe beam (i.e. in the direction of the pump beam) to that of the incoming 
probe beam. From eqs. fl63l) and fl65l) . this is: 

r^= |X(-i)|2 = j2(5). (72) 

For ease of presentation of our results, it is also convenient to define quantities rj' and g'. 
These quantities are respectively analogous to r] and g, but with the roles of the pump and 
probe beams exchanged. 



V. RESULTS 



A. Analytical study of the behavior of intensities. 

First we discuss diffraction from a single beam. The transmitted energy is divided between 
beams of different orders n. The amplitude of diffracted beams is given by |X*^"^p oc Jn{B) 
(eq. (l63|) ). where we recall, from eq. (l58!) that B is the amplitude of the additional phase 
variations induced by the spatial modulations in the liquid crystal layer. We introduce the 
energy conservation parameter 

k=n 

Sn{B) = JliB). (73) 

k=—n 
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The quantity Sn{B) describes the proportion of transmitted energy that is distributed be- 
tween diffracted beams of orders from — n to n. We note that in the hmit n ^ oo necessarily 

Sn 1- 

In FiglHwe show the dependence of on the phase modulation parameter B for low n. 
For small B {B < 1), Si{B) ^ 1. In this regime almost all transmitted energy is either in 
the directly transmitted beam, or in the first-order diffracted beam. This is the regime in 
which energy exchange between two beams is most effective. 

As B increases, an increasing proportion of the transmitted energy is transfered to outly- 
ing diffracted beams. The quantity S2 remains essentially unity until B ^ 1.5, while S3 only 
noticeably departs from unity at i? ~ 2.5. We shall return to the problem of the asymptotic 
behavior of for large n,B elsewhere. For this study, however, we shall be interested 

in the small B regime. 




B 

FIG. 4: The dependence of energy conservation parameter on phase modulation parameter 

B. Solid line: si, dashed line S2 and dotted line S3. See text for further discussion. 



We now analyze the effect of energy exchange in the presence of both incident beams. We 
shall take equal intensities in eq. fITU]) for incident beams m = 1. In principle the parameter 
m [eq. fl69p ]. measuring the ratio of the intensities of the beams, can take any value. In our 
calculations we shall suppose m = 1; this corresponds to equal intensity beams. This will 
enable us to make contact with previous studies which have also used this value. In 
addition it is easy to monitor energy transfer between beams. We note that in devices we 
may well expect that m -C 1, so that a large reservoir of pump beam energy is available to 
amplify a given probe beam. 
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The degree of energy transfer is critically dependent on the quantity A, defined in eg. (1581) . 
When A = 0, the modulation of the phase retardation is in-phase with the intensity mod- 
ulation due to the beam interference. If m = 1, eq.l lTOl) implies that there is no net energy 



exchange between the beams. This is consistent with general intuition 5 



that a phase 



difference between the intensity and dielectric modulations is required for beam coupling. 
Interestingly, although we do not pursue this here, this rule no longer holds for the m 7^ 1 
case. 

The maximum energy transfer between beams occurs when sin(A) = 1. In this case the 
two outgoing beams obey the following rule: 

Ji oc {MB)-Ji{B)f 

h oc (Jo(5) + Jl(5))^ (74) 

The behavior of these functions is shown in FiglSl The function hiB) (solid curve on the 
graph) reaches a maximum value of ~ 1.48 at B ^ 0.85. We also plot Ii{B), and note that 
this reaches a minimum(at zero) for B ^ 1.4. The quantity 

m=^-ihiB) + hiB)) (75) 

denotes that proportion of the energy of the incident beams which remains in the two initial 
beams directions. The quantity I{B) is unity for 5 = (at which there is, however, no 
energy exchange) and reduces steadily with B. Close to the maximum l2{B) ^ 1.48 at 
B ^ 0.85, I{B) ^ 0.8, reducing monotonically to I{B ^ 2) ^ 0.4. However, for B < 1, 
Si{B) is essentially unity. The energy lost from the primary beams reappears as other |n| = 1 
diffracted beams. Subsequent maxima of hiB) take values less than unity, in regimes in 
which a substantial proportion of the transmitted energy is lost in outlying diffracted beams. 
Thus the appearance of energy transfer between beams is restricted to the first maximum 
ofhiB). 



B. Dependence on external parameters 

We now examine quantitatively the energy exchange process, using experimentally plau- 
sible parameters. A list of parameters in the problem is given in Table [H together with 
meanings of these parameters, and where appropriate, the numerical values that we have 
used in our calculations. 
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B 

FIG. 5: The behavior of functions /2(-B)(solid curve) and /i(B)(dot curve) when sin(A) = 1. 

TABLE I: Table of parameters 



Parameter 


Value 


Description 


A 


0.63/im 


Wavelength of incident beams 


L 


20fim 


Thickness of the film 


£\\ 


1.52, 1.72 


Dielectric permittivities of the liquid crystal 




variable 


Angle of propagation inside liquid crystal 


6 


variable 


Phase shift between the interference patterns at top and bottom surfaces. 






This is the surrogate for the angle of incidence which we do not include 






explicitly. 


7 


variable 


Half-angle between beams defining the dimensionless grating wave- 






vector fi = qL = 2A;Ly/|^cosV'sin7 




variable 


Non-dimensional grating wave vector /x = qL. For 7 ~ 2.4*^, ~ 6. 


V 


1 


Dimensionless voltage = LEq{-^)^''^ 

K 


m 


1 


Ratio of intensities of incoming beams 



The final element in the theory enabling comparison with experiment is the response 
of the surface potential $1 to the local beam intensity. We do not however have a micro- 
scopic photoelectrochemical theory to describe this process. In principle $1 is a measurable 
quantity, although in practice the measurement may be difficult to carry out. 

In our initial calculations we suppose $i/$o = 1 (eq j271) and the external voltage u = 1. 




24 



We first investigate tfie dependence of energy exchange effect on We set the grating 
period = 6, which corresponds to a grating wavelength A equal to the cell thickness L. 
To see the influence of both the in-phase and out-of phase components we choose the phase 
shift of the interference patterns 6 = 7r/2. 




0.5 1 

■tjj, Rad 



FIG. 6: Functional dependence of the phase modulation parameter B as a function of the internal 
angle ip. Fixed parameters: 5 = tt/2, fi = 6 and <I>i/<^o = 1- 

The key intermediate parameters governing g [eq.l lTTi) ] and r] [eq. (172l) ] are the phase 
modulation parameters B and A. In Fig. [6] we plot the phase modulation parameter B 
as a function of ip, the angle of propagation inside the hquid crystal. The principal result 
from Fig. [6] is that B < 1 everywhere, which implies that only transmitted and first order 
diffracted beams occur. In Fig. [7|we plot It can be seen from eq. (170|) that the phase 

modulation parameter A governs the sign of the energy exchange. For < A < vr sin A > 0, 
and the probe beam is amplified by the pump beam. However for tt < A < 27r sin A < and 
the pump beam is amplified. In Fig. [6]i?(^) = for ip ^ 0.68. At this point the quantities 
A [eq. fl59|) ] and C [eq. fl60l) ] are both equal to zero and change sign. This implies a sudden 
phase shift of vr in the phase modulation parameter A, which indeed occurs at ip ^ 0.68 in 
Fig. El 

The behavior of the gain is shown in Fig. [HI Maximal gain is achieved when smA{ip) = 1, 
corresponding to ip ^ 0.95 (see eqs.ljTO j [741) ). At ip = the gain g{ip) = 1. We also note that 
at this point, from the symmetry of the system there is no energy exchange. The quantity 
/ plotted on this graph is almost unity everywhere. This means that there are no energy 
losses and all incident energy is distributed between outgoing probe beam and outgoing 
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FIG. 7: Functional dependence of the phase modulation parameter A as a function of the internal 
angle. Other parameters are fixed: 5 = tt/2, fi = 6 and ^i/^q = 1. 

pump beam. 

2 
1.5 

9 

1 

0.5 



"0 0.5 1 

ip, Rad 

FIG. 8: Energy exchange between beams as a function of the internal angle ip. Solid line: gain of 
probe beam. Dashed line: gain of pump beam. Dotted line: / (see eq. (j75p ) represents the degree 
of energy conservation in the system. Fixed parameters: 5 = tt/2, fi = 6 and $i/<^o = 1- 

In Fig. M we plot the pump beam diffraction efficiency. This measures the proportion 
of energy diffracted from the pump beam in the probe beam direction. The diffraction 
efficiency 77 is a function only of B [eq. fl72|) ]. The maximum possible diffraction efficiency is 
given by rji^B) ^ 0.338 at i? ~ 1.83. This maximum value does not depend on the details 
of our model and remains true for thin gratings [l^. But in Fig. [9], the maximum value of 
diffraction efficiency t] ~ 0.08 occurs at ip ^ 1.05. The discrepancy between the theoretical 

26 




0.08 - 



0.04 




FIG. 9: The dependence of diffraction efficiency rj on the internal angle ip. Fixed parameters: 
S = 7r/2, ;u = 6 and $i/$o = 1- 

maximum and this maximum can be ascribed to the fact that B{iIj) < 1 everywhere. 

The phase shift 6 [eq.Q] also strongly affects the energy exchange characteristics. In 
Figs. [MH 6 = it/2 and the modulation of the dielectric and of the energy along the cell 
are out-of-phase. In Fig. [10], we put 6 = tt; now these modulations are in-phase. Now it 
is possible to transfer energy from the probe beam to the pump beam, by contrast with 
previous case in which negative energy transfer never occurs. Thus, whereas in Fig. [H]^' > 1 
everywhere, in Fig. [TOl g — 1 can take either positive or negative signs. 

We now turn to the study of energy exchange as a function of the angle between interfering 
beams. Increasing the angle 27 between beams increases the non-dimensional wave vector 
fi = qL. Quantitatively, choosing parameters given in Table [U we find that for 7 ^ 2A^, 
fi ^ 6. We examine the dependence g{fi) for tp = 0.95, i.e. at the maximum of g appropriate 
to Fig. [HI The dependence of the gain on fi is shown in Fig. [TTl For fi < A, B is larger 
than unity (see inset a)) and J < 1 (see eg. (1751) ) is less than unity. In this regime there is 
considerable energy loss due to probe beam diffraction into diffraction orders of order n > 1. 
However for fi > 4, B < 1 and the quantity / is close to one. In this regime energy is 
conserved. The maximal gain is achieved at /i ~ 6. This maximum is consistent with our 
results from Figs. [MH 

In Fig. [12] we plot the diffraction efficiency rj. The maxima of 0.338 occur at /i ~ 0.65 
and fi ~ 2.65, corresponding to maxima in the Raman-Nath regime 16|]. But in our case 
these maxima occur for /i < 4, where as we have seen above, there is considerable diffractive 
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FIG. 10: As for Fig. [H but with 6 = tt (see text). Energy exchange between beams as a function 
of tp. Sohd hne: gain of probe beam. Dashed hne: gain of pump beam. Dotted hne: / (see Fig. 
[8|). Insets: a) Phase modulation amphtude S as a function of the mean angle of propagation ■0; 
b) Phase modulation parameter A{ip). 
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M 

FIG. 11: Energy exchange between beams as a function of /x, the non-dimensional grating wave 
vector. Solid line: gain of probe beam. Dashed line: gain of pump beam. Dotted line: / (see 
Figs. [HI nop . Insets: a) Phase modulation amplitude B as a function of fi; b) Phase modulation 
parameter A{p). Fixed parameters: 6 = 7r/2, ip = 0.95 and <I>i/<I>o = 1- 

energy loss. The physical relevant maximum occurs for T]{fi = 6.5) ~ 0.07. Here there is 
insignificant energy loss. 

$1 

The strength of the grating modulation (eg. (1271) ) depends on the ratio — . In all our 

<±>o 

$1 

previous plots (Figs. [6lfT2l) we have set the ratio — = 1. This ratio can be modified in two 

$0 
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FIG. 12: The dependence of diffraction efficiency r] on the non-dimensional wave vector /x. Fixed 
parameters: 6 = vr/2, ip = 0.95 and <I>i/<I>o = 1- 

ways. In principle one can change $i by changing the surface preparation. Alternatively 
(and more simply) one can apply an external voltage $o across the liquid crystal. Here we 
suppose $1 and $o to be independent quantities. 

We now investigate the effect of external voltage on the energy exchange for g^ip = 
0.95, yU = 6). This corresponds to the point (see FigJS]) where g is maximal with respect to 
varying ip with other parameters as taken in Table [B The optical modulation is a strong 
function of the director modulation. It is thus useful to analyze the director modulation as 
a function of external applied field. 

We first make a quantitative analysis of the behavior of the director modulation as a 
function of voltage. The liquid crystal director distribution is given by eq.( !271) : 



The voltage enters this expression explicitly through the multiplier — , and implicitly 

L 

through the quantity k, where = /i^ + i/^, z/ = — = $o(~")"'^^^ is a rescaled voltage. 

In the limit of high voltages, the n term in eq. (I27p can be neglected, and hence the 
reorientation 9{x,z) oc $o^. Hence at high voltages the director approaches a uniform dis- 
tribution, in which case the beam-coupling disappears and g = l- In our approximation, the 
low voltage limit $ ^ is inaccessible; the minimum voltage for which our approximations 




(27) 



$1 



29 



1.6 
1.4 
1.2 

9 1 
0.8 
0.6 
0.4 
0.2 


FIG. 13: Dependence of gain of probe beam as a function of <&o/*^i- This is equivalent to increasing 
the external field. Note that the abscissa starts at 1, the lowest value for which our treatment is 
valid, although the difference between 1 and is not visible on this scale. Further structure may 
occur in the region ^o'^'i ^ 1- 

are valid is $o ~ 1- This follows because we have assumed that the modulated component of 
the electric field is small by comparison with the external field in the derivation of the liquid 
crystal director distribution eq. (!24l) . The dependence of gain on the external potential is 
shown in Fig. [131 For $o = 1 fl' ~ 1-3. The gain reaches a maximum g 1.44 at $0 ~ 8- For 
$0 > 8 the energy exchange parameter g decreases monotonically toward a value of unity 
(i.e. no energy exchange) in the high potential limit. 

VI. DISCUSSION AND CONCLUSIONS 

In this paper we have carried out a model phenomenological calculation of energy ex- 
change between beams incident on a thin liquid crystal grating sandwiched between two pho- 
toconducting layers. In this model calculation, the liquid crystal is subject to homeotropic 
boundary conditions, but this is not an essential feature of the model. The energy ex- 
change involves diffraction by an induced grating, with the exchange occurring when each 
incident wave is diffracted into the outgoing path of the other. We find that there is a 
regime in which significant energy exchange can occur, without leakage into other higher 
order diffracted waves. There is also another regime in which such leakage does occur. 

We find a maximal gain oi g = 1.45, and this occurs for a grating wavelength of the order 
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of the thickness of the sample. This appears to be a robust result, and qualitatively consistent 



with experiment 



26l ] . There is also significant dependence on the bulk voltage. In the limit 



of high voltage there is no effect (because there is no director modulation and hence no 
grating). As the voltage is reduced the effect increases. Our calculation does not permit the 
evaluation of the low- voltage limit, but we do find a maximum when the ratio of the external 
voltage to the surface modulation is about 10. Unfortunately the phenomenological nature 
of our calculation does not permit us to make any quantitative predictions with respect to 
actual voltage or beam intensities required to achieve this. However, the gain maximum as a 
function of voltage is also manifested as a maximum with respect to varying beam intensity. 
Although this calculation is vague with respect to quantitative prediction, we believe that 
the existence of a maximum as a function of voltage is a qualitatively robust result. 

The calculation is broken down into a number of parts. Firstly we have supposed that 
interference between the incident beams affects the photoconducting layers by only changing 
the electric potential at the boundaries of the sample, and does so in proportion to some 
power of the beam intensities. Secondly we have calculated the modification of the electric 
field inside the liquid crystal sample, supposing that there is a zeroth order field due to some 
imposed bulk potential. Thirdly, we have used the electric field to calculate the modulated 
director distribution, which necessarily then acts as an optical grating. Fourthly, we have 
investigated the transmission of each beam through the modulated liquid crystal layer using 
a WKB-like approximation in the spirit of geometrical optics. The result of this calculation is 
a phase and amplitude optical profile for the extraordinary wave along the outgoing surface. 
Finally, using this surface optical profile we have used the Kirchoff method to evaluate 
the far field, and hence diffraction and inter-beam energy exchange. We may note that the 
resulting expression for the intensity of the diffracted beams of different orders is reminiscent 
of the analogous calculation for diffraction through a thin grating in the Raman-Nath regime 
[igI i . However, the standard Raman-Nath calculation does not apply here. The existence of 
inhomogeneities in the dielectric function perpendicular to the grating direction complicates 
the calculation. 

Some features of our calculation are simplified in order to make the problem tractable. We 
have measured the potential induced at the surface with respect to the bulk voltage across 
the cell. An alternative low external voltage expansion would also in principle be possible, 
but we have not pursued this approach here. In addition, we have solved the liquid crystal 
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director in a one-elastic-constant approximation, and also linearized the Euler-Lagrange 
equations coupling the director and the electric potential. Neither of these approximations 
is essential, but there is significant potential advantages in obtaining analytic formulas in 
what could otherwise be a computational minefield. The optical scattering in the system 
itself implies that the picture in which the incident beams penetrate the liquid crystal layer 
unimpeded to provide potential modulations at the outgoing surface is only the first step in 
an iterative procedure. While it is possible to carry out this iteration in our model, we have 
chosen not to do so. This is partly because we would lose what analytic simplification we 
have achieved. Also, however, given that we only have a phenomenological model, we would 
in any case be no closer at this stage to a quantitative comparison with experiment. 

One particularly interesting feature of our calculations is that we are able to identify 
separately components of the refractive index modulation which are respectively in-phase 
and out-of phase with the mean optical field intensity. It is often stated that if the re- 
fractive index and optical field modulations are in-phase with respect to each other, then 
no two-beam coupling would be expected. However, notwithstanding the inaccuracies and 
approximations involved in our calculations, we find that this statement is unambiguously 
false. The in-phase beam coupling is indeed lower, but no by means identically zero. 

Although we have been able to obtain a semi-analytic form for the beam coupling, the 
calculation is complicated. It involves electric fields, director distributions and light trans- 
mission through an inhomogeneous medium. The result is that the final magnitude of the 
effect under consideration seems to bear no simple relation to the rather large number of 
parameters which enter the problem. We can say definitively that the external voltage, the 
angle of incidence, the angle between the two beams, not to mention the thickness of the 
cell and the liquid crystal elastic constants, all play an important role, and furthermore 
the response is not monotonic. From an engineering point of view, there are clearly several 
possible ways to control the energy exchange process. 

But apart from the pronounced maximum in beam coupling when the grating width is of 
the order of the thickness of the sample, we are unable at this stage to make further robust 
comments concerning the functional relationships without resorting to specific calculations. 
We cannot say whether further studies, and in particular a reliable microscopic theory, will 
clarify the situation. 

The theory presented in this paper can be developed in a number of ways. It can be 
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trivially extended to liquid crystal cells in which only a single photoconductive layer at- 
tached. Alternatively, we might extend the present work, involving homeotropic surfaces, 
to low voltages, or to liquid crystal cells with homogeneous boundary conditions. Such a 
theory would be applicable, for example, to the experiments of Pagliusi and Ciparrone 
We also note that in this paper, the director distribution throughout the sample is clustered 
around the homeotropic direction. However, one might expect intuitively that the most 
dramatic effects would occur when the voltage modulation and the external field conspire to 
produce large director shifts between one part of the sample and another. Our framework 
may permit such a calculation. 

The main weakness of the theory concerns the nature of the relationship between the 
potential modulations $i and the beam intensities. The lack of relevant experimental data 
is partly because, as far as we are aware, the present paper is the first suggestion that the 
main mechanism for photorefractive beam coupling involves this quantity. We are hopeful 
that future work will therefore remedy this deficiency. An experiment which measures this 
quantity might involve Frederiks transition measurements in the presence of an externally 
applied optical beam. 

At a later stage, we would also hope to make contact between this theory and a more 
microscopic theory which elucidate processes in the photoconductive media, and at the 
photoconductive layer-liquid crystal interface. A second weakness involves the geometric 
optics approximation, and this restricts our calculations to the short wavelength limit. In 
most liquid crystal cells, this will be sufficient, but in principle longer- wavelength corrections 
are interesting. Indeed, we are currently carrying out optical calculations in which the optics 
is treated by solving the Maxwell equations exactly. Such a scheme will automatically permit 
a self-consistent solution of the optics-potential-elastic problem. 
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